To Appear in ApJ 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



KEPLER EXOPLANET CANDIDATE HOST STARS ARE PREFERENTIALLY METAL RICH 

Kevin C. Schlaufman 1 and Gregory Laughlin 

Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064 
Received 2011 February 17; accepted 2011 June 23 

ABSTRACT 

We find that Kepler exoplanet candidate (EC) host stars are preferentially metal-rich, including the 
low-mass stellar hosts of small-radius ECs. The last observation confirms a tentative hint that there is 
a correlation between the metallicity of low-mass stars and the presence of low-mass and small-radius 
exoplanets. In particular, we compare the J — H — g — r color-color distribution of Kepler EC host 
stars with a control sample of dwarf stars selected from the ~150, 000 stars observed during Ql and 
Q2 of the Kepler mission but with no detected planets. We find that at J — H = 0.30 characteristic 
of solar-type stars, the average g — r color of stars that host giant ECs is 4-er redder than the average 
color of the stars in the control sample. At the same time, the average g — r color of solar-type 
stars that host small-radius ECs is indistinguishable from the average color of the stars in the control 
sample. In addition, we find that at J — H = 0.62 indicative of late K dwarfs, the average g — r color 
of stars that host small-radius ECs is 4-er redder than the average color of the stars in the control 
sample. These offsets are unlikely to be caused by differential reddening, age differences between the 
two populations, or the presence of giant stars in the control sample. Stellar models suggest that 
the first color offset is due to a 0.2 dex enhancement in [Fc/H] of the giant EC host population at 
ill* « 1 Mq, while Sloan photometry of M 67 and NGC 6791 suggests that the second color offset is 
due to a similar [Fe/H] enhancement of the small-radius EC host population at ps 0.7 Mq. These 
correlations are a natural consequence of the core- accretion model of planet formation. 
Subject headings: planets and satellites: detection — planets and satellites: formation — protoplane- 
tary disks — stars: abundances 



1. INTRODUCTION 



The groundbre aking early 
Kepler mission ( Borucki et al 
2010a l Uenk ns et all 12010a ; 



Brvson et all l201fl iBatalha et al.l l2010al lbl: lHaas et all 



o f NASA's 
IKoch et all 
Caldwe 1 et al.l 120101: 



r esults 
[20105 



2010t IJenkins et all l2010bl) have spectacularly added 
to our understanding of the frequency and prope r- 
ties of planets i n the Galaxy jlBorucki et all [201 lallbl: 
Ford et all 120111: iLissauer et all 1201 lbt iMoorhead et al l 
20111: iHoward et all 120111 : lLatham et all 120111 ). Eight 
planetary systems discovered by Kepler have already 
been confirmed by radial velocity or transit timing: 
four hot Jupiters (IKoch et ail l2010bl: iDunham et all 
120101: lLatham et ail 120101 : IJenkins et all l2010cD : two 
Saturn- mass planets and a su per-Earth in the Kepler-9 
system dHolman et all l2010f): a Neptune-mass planet 
Kepler-4b |Borucki et alJl2010bft: Kepler- 10b the small- 
est exoplanet yet known (|Batalha et all 1201 ID ; and the 
closely-aligned six tr ansiting planet system Kepler- 11 
(|Lissauer et al.ll2011a| ). 

The properties of the exoplanet candidate (EC) pop- 
ulation discovered by Kepler will provide strong con- 
straints on models of planet formation, especially for 
low-mass stars and small planets. It is well established 
that metal-rich solar-type stars are more likely to host gi- 
ant planets than more metal-poor solar-typ e stars (e.g., 
iSantos et all l2004t iFischer fe Valentil |2005[ ). Likewise, 
metal-rich low-mas s stars also seem more likely to host 
giant planets (e.g.. Uohnson fc A pps 2009). Meanwhile, 
preliminary results indicate that near solar metallicity 



there is not much of a correlation between solar-type host 
star metallicity an d the likelihood o f hosting Neptune- 
mass planets (e.g. lUdrv et all 120061: iSousa et all 120081: 
iBouchv et a"ill2009D . There are also trends in host stellar 
mass, as higher mass stars are mor e likely to host giant 
plane ts than lower mass stars (e.g.. Uohnson et"aT1l2007L 
l2010l: iBorucki et alJl2011bH 

In [Schlauf man fc Laughlinl (|2010D . we tentatively 
noted a hint that M dwarfs that host low-mass plan- 
ets are more likely to be metal rich than expected based 
on random sampling from the field M dwarf population. 
We explained this correlation in the conte xt of the core- 
accre t ion model of planet formation (e.g.. iPollack et al l 
T996] : llda fc Linl 12001 iLaughlin et all 120041: llda fc Lir 



2005tlHubickvi et al.ll2005l : iMordasini et al.ll2009allbH In 
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particular, we argued that a more metal-rich protoplane- 
tary disk will almost certainly have a higher surface den- 
sity of solids. That increased solid surface density may 
enable the rapid formation of the several Earth-mass em- 
bryo necessary to form the core of an ice giant and ac- 
crete gas from the protoplanetary disk before the gaseous 
disk is dissipated. That correlation might also indicate a 
lower limit on the amount of solid material necessary to 
form planets. 

In this paper, we determine whether Kepler EC host 
stars are more metal rich than stars with no detected 
planets. To that end, we use a metallicity-sensitive color- 
color plot to examine the optical and infrared properties 
of the sample of Kepler EC host stars relative to a control 
sample of dwarf stars observed by Kepler but that were 
not observed to host transiting planets. We describe our 
analysis procedures in Section 2, we discuss the results 
and implications of our analysis in Section 3, and we 
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summarize our findings in Section 4. 

2. ANALYSIS 

We are interested in using the available photometric 
data for stars in the Kepler field to determine if there is 
an offset in metallicity between the Kepler EC host popu- 
lation and the population of stars observed by Kepler but 
with no detected planets. In particular, we compare the 
mean g — r color of Kepler EC hosts with the mean g — r 
color of stars observed by Keple r but with no dete cted 
planets at constant J — H color. Ilvezic et al.l (|2008l ) ex- 
plored the potential of SDSS photometry to produce ac- 
curate photometric metallicities for metal-poor stars in 
the halo of the Milky Way. Unfortunately, the accuracy 
of those methods depends crucially on the availability 
of precise SDSS M-band photometry. Nevertheless, even 
when u-band data is lacking, <?,r, and i photometry can 
still provide ver y useful photometric metallicities (e.g., 
I An et al.|[2009bT) . 

Our problem is made easier because we are only in- 
terested in the relative (not absolute) metall icities of the 
two p opul ations. At the sam e time, unlike Ilvezic et al.l 
(2008) and I An et all (|2009bD . we have the added benefit 
of 2MASS JHK photometry in the Kepler field and ap- 
parent magnitude range. As a result, we can use infrared 
colors (e.g., J — H) as a proxy for effectiv e temperature 
r eff . The Kepler Input Catalog (KIC - iBatalha et all 
l2010bt iBrown et alj I2011D provides reasonably accurate 
estimates of log g, so we can ensure that giant stars do not 
contaminate our control sample of stars with no detected 
planets. Stars with M* w 0.7 Mq have main sequence 
lifetimes much longer than 10 Gyr, and because age has 
little effect on the colors of low-mass stars on the main 
sequence, age differences between the two populations (if 
they exist) are unlikely to affect our analysis. Since we 
compare dwarf stars with similar T c ff and logg, the only 
stellar parameter that can vary between the Kepler EC 
hosts and the control sample is [Fe/H]. As metallicity is 
the only stellar parameter that can vary, the best expla- 
nation for systematically red g — r colors in the EC host 
population relative to the control sample is that the EC 
host population is enriched in metals. 

2.1. Color Offsets 

For each of the 997 stellar hosts of the 1,235 Kepler ex - 
oplanct candidates announced in lBorucki et al.l (|2011bf ). 
we obtain all available broadband photometr y (includ- 
ing 2MASS JHK from ISkrutskie et all 12051 and red- 
dening data from the KIC. Simultaneously, we obtain 
equivalent photometry and reddening data from the same 
sources for a random sample of 10,000 dwarf stars from 
the ^150,000 stars observed during Ql and Q2 of the 
Kepler mission with no candidate planets and KIC-based 
log g > 4. This control sample is subject to the same se- 
lection effects that were applied in the Kepler field to pro- 
duce a list of stars to search for transiting planets. Con- 
sequently, differences in the characteristics of the stars in 
the control sample and the stars that host ECs are not 
related to the selection effects applied in the Kepler field 
to identify a sample of stars suitable for transit observa- 
tions. 

Metal-rich stars typically have redder optical colors 
than solar metallicity stars, as the forest of iron lines 



in the atmospheres of metal-rich stars preferentially ab- 
sorb blue photons. At the same time, there are other 
factors that affect the color of a star, most notably red- 
dening due to the interstellar medium, its age, and its 
evolutionary state (e.g., dwarf or giant). 

In Figure [TJ we plot a J — H — g — r color-color plot 
for the control sample and the EC host sample. We sub- 
divide the EC host sample into two sub-samples: those 
stars that host at least one giant EC (e.g., R p > 5 i?®) 
and those s tars that host no g i ant ECs. We use P adova 
isochrones ([Marigo et al.ll2008t iGirardi et alJ[2010l ) for a 
2 Gyr stellar population at [Fe/H] = -0.1 and 0.1 to de- 
termine the effect of metallicity on the stellar locus in a 
J—H — g — r color-color plot. The yellow arrows show the 
effect of increasing [Fe/H] by 0.2 dex on g — r at constant 
J — H. Note though that the yellow arrows do not con- 
nect stars of constant mass, as a 0.2 dex [Fe/H] enriched 
star will be about 5% more massive than a solar metal- 
licity star at constant J — H color. We do not use the 
Padova isochrone redward of J — H = 0.52, as the theo- 
retical stellar models and atmosphere s become unreliable 
after that point (e.g,. IAn et al. 2009a) . Conseq u ently, we 
use the M 67 fiducial sequence from lAn et al.l (|2008j ) to 
illustrate the morphology of a solar metallicity popula- 
tion at red J — H color. Age can also affect the color 
of main sequence star, though near J — H = 0.62 age 
has little effect on the main sequence. Giant stars have 
slightly bluer g — r colors than dwarf stars at constant 
J — H, so the presence of giant stars in the control sam- 
ple could in principle produce a color offset between the 
Kepler EC host population and the control sample. For- 
tunately, KIC-based log g estimates are accurate enough 
to ensure that g iant stars are at most a few percent of our 
control sample (jBasri et alJl2011b [Brown et al.ll2011l ). 

Two trends are immediately apparent in Figure [TJ 
First, at J — H > 0.6 typical of late K dwarfs, giant 
ECs become very rare relative to smaller ECs. We in- 
terpret this as evidence in support of the correlation be- 
tween host stel lar mass and likelihood of hosting a gi- 
ant pla net (e.g.. Uohnson et aX1l2007l 120101 iBorucki et al.1 
I2011b| ). Second, at J - H » 0.62 dominated by late K 
dwarfs, ECs are preferentially found around the reddest 
part of the stellar locus as defined by the control sam- 
ple. We quantify the significance of this latter trend in 
Figure [2j 

To determine the significance of the observation that 
low-mass stars that host ECs are on average redder in 
g — r than the stellar locus as defined by the control 
sample, we bin both the control sample and each sub- 
sample of EC hosts into 0.16 mag bins in J — H color 
and compute the mean and median g — r color in that 
bin. We estimate the error in each mean and median 
by bootstrap resampling. We find that at J — H w 0.3 
characteristic of solar-type stars, giant EC hosts have 
preferentially red g — r colors at the 4-cr level relative 
to the control sample. On the other hand, at that same 
J — H color, the hosts of small ECs are indistinguishable 
from the control sample. In addition, at J — H « 0.62 
the stars that host small ECs are 4-cr redder in g — r 
than the control sample. We summarize the results of 
our calculations in Table [TJ 

2.2. Possible Reasons for Color Offsets 
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We argue that the significant color offsets apparent in 
Figure [2] and Table [1] are the result of the enriched metal- 
licity of the EC host sample. Indeed, it reassuring that 
we identify the hosts of giant ECs as significantly red in 
g — r at constant J — H, as it is well established that the 
stella r hosts of giant planets are preferentially m etal rich 
(e.g.. lSantos"etaL|[200l [Fischer fc Valentill2005[ ). Other 
possible explanations for the observed offset are selection 
effects, differential reddening, systematic age differences 
between EC hosts and the control sample, and the pres- 
ence of giant stars in the control sample but not the EC 
host sample. 

2.2.1. Selection Effects 

The depth of a transit is proportional to (i? p /i?*) 2 , 
so it is easier to identify the signal of transiting planet 
around a small star. On the main sequence, both age 
and metallicity can affect the radius of a star. Accord- 
ing to the Padova isochrones for a 2 Gyr population 
(jMarigo et al.ll200l iGirardi et aLll20Toh . at J - H = 0.3 
a star with [Fe/H] = 0.1 is about 5% larger than a star 
with [Fc/H] = —0.1 and the same J — H color. At 
J-H = 0.62, a star with [Fe/H] = 0.1 is about 7% larger 
than a star with [Fc/H] = —0.1 and the same J — H color. 
According to the Padova isochrones, at J — H = 0.3 a 
5 Gyr old star is about 6% larger than a 2 Gyr old star 
at the same J-H color. At J - H = 0.62, a 10 Gyr 
old star is about 2% larger than a 2 Gyr old star at the 
same J — H color. Consequently, if there were no correla- 
tion between host metallicity and probability of hosting 
a planet or possible tidal destruction of aged exoplanet 
systems, transit surveys would be more likely to identify 
transiting planets around young, metal-poor stars. As 
we show in Section 2.2.3, age on the main sequence has 
very little effect on the g — r color of a star at constant 
J — H, so these effects are not likely to produce the color 
offsets we observe. 

The g — r color offset at constant J — H color that we 
observe between the Kepler EC host sample and the con- 
trol sample is only meaningful if the two samples were 
subject to the same selection effects. As we argued in 
Section 2.1, both samples were selected according to the 
same Kepler target selection algorithm. To verify that 
the two samples are similar, we plot in Figure [3] six color- 
color plots and tw o color-magnitude diagrams (similar to 
ICovev et al.ll2007h . In all cases, the control sample out- 
lines the distribution of both the sample of giant EC 
hosts and the sample of small-radius EC hosts. The EC 
hosts that are outliers in each panel may be unequal mass 
binary star systems or systems with poorly-estimated 
reddening. The Kepler IDs (KOI Numbers) of a few 
of our most extreme color outliers are: 1161345 (984), 
2446113 (379), 5356593 (644), 8162789 (521), 10470206 
(335), 10514430 (263), and 11465813 (771). 

In short, the age or metallicity of the host star of an 
exoplanet does not significantly affect the probability of 
detecting a transiting planet, and our control sample is 
a reasonable control sample for our measurement of rel- 
ative color offsets. We therefore conclude that selection 
effects are unlikely to produce the color offsets we observe 
in Figure [2] and Table [U 

2.2.2. Differential Reddening 



Reddening due to the interstellar medium can affect 
the observed color of a star. In our differential analy- 
sis, this will only affect our result if the reddening of the 
control sample is systematically different than the red- 
dening of the Kepler EC host sample. Fortunately, Fig- 
ure |4] indicates that control sample is subject to the same 
reddening distribution as the EC host sample. Like the 
EC host sample, the control sample is also spread more 
or less uniformly over the Kepler field, so any angular 
dependence on reddening is unlikely to produce a signifi- 
cant color offset between the control sample and the EC 
host sample. We therefore conclude that differential red- 
dening between the control sample and the Kepler EC 
host sample is unlikely to produce the color offsets we 
observe in Figure [2] and Table [TJ 

2.2.3. Systematic Age Differences 

As stars age on the main sequence, their colors can 
evolve. According to the Padova isochrones, at J — H = 
0.3 a 5 Gyr old star is about 0.1% redder in g — r than a 
2 Gyr old star at the same J — H color. At J — H = 0.62, 
a 10 Gyr old star is about 0.5% bluer in g — r than a 2 
Gyr old star at the same J — H color. Clearly, colors of 
stars do not evolve significantly on the main sequence. 

Tidal evolutio n in exoplanet systems can destroy close- 
in planets (e.g., iGu et all [20031 iMardling fc LinT[200l . 
Transit surveys are strongly biased towards close-in sys- 
tems, so it is possible that transit surveys may be biased 
towards young stars, as tidal evolution may have already 
destroyed planets that once orbited older stars. Assum- 
ing that planets arrive in the close proximity of their host 
star on nearly circular orbits and that the timescale for 
eccentricity damping is very short, the n the timescale for 
tidal disruption is approxima tely (e.g., Ilbgui fc Burrows! 
[2001 ISchlaufman et al.l[20T?l 

r dis - — Gl/2 ^ — [l-(R*/a ) /j, (1) 

where ao is the initial semimajor axis of the planet before 
tidal evolution, Q' t is the specific dissipation function of 
the host star, G is Newton's gravitational constant, M p 
is the mass of the planet, and R* is the radius of the 
host star. If tidal evolution efficiently destroys planets 
in the Kepler sample, then the hosts of Kepler ECs may 
be preferentially younger than stars in the control sam- 
ple. The median period of giant ECs is 13 days, while 
the median period of small-radius ECs is 10 days. As we 
plot in FigureEJ the timescale for tidal disruption of such 
systems is in excess of 100 Gyr. At the same time, as we 
showed above, age has little effect on the g — r color of a 
main sequence star at constant J — H. We therefore con- 
clude that systematic age differences between the control 
sample and the Kepler EC host sample or tidal evolution 
in EC systems are unlikely to produce the color offsets 
we observe in Figure [5] and Table [TJ 

2.2.4. Presence of Giant Stars 

Giant stars have stellar radii i?„ ~ 100 i?o, and be- 
cause transit depth is proportional to (R p / R*) , the tran- 
sit of a planet with radius R p in front of a giant star with 
i?* ~ 100 Rq is about 10,000 times harder to detect than 
the transit of the same planet in front of a dwarf star with 
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radius J2* ~ 1 Rq. As a result, Kepler is unlikely to iden- 
tify transiting planets around giant stars, so the sample 
of Kepler EC hosts is very likely free of giant stars. 

While we only used stars with KIC-based log g > 4 in 
our control sample, it is possible that a few p ercent of the 
stars in our control sam ple are giant stars ([Basri et al.l 
l201lHBrown et al.ll201ll) . Giant stars have slightly bluer 
g — r colors at constant J — H color than dwarf stars, 
so the presence of a significant number of giant stars in 
our control sample could explain the 0.08 mag g — r color 
offset we observe at J — H — 0.62 between the sample of 
Kepler EC hosts and our control sample. 

To investigate the affect of significant giant star con- 
tamination in our control sample, we obtain all available 
broadband photomet ry (including 2MASS JHK from 
ISkrutskie et a,l.|[2006)f) and reddening data from the KIC 
for a sample of stars observed by Kepler during Ql and 
Q2 with KIC-based log g < 3. The stars in this sam- 
ple are likely giant stars. Recall that the 10,000 stars 
in our control sample have KIC-based logg > 4, so we 
can add stars from the sample of likely giant stars to 
our original control sample of likely dwarf stars to create 
giant star contaminated versions of our our control sam- 
ple. In particular, we create four versions of the giant 
star contaminated control sample with differing levels of 
contamination: 30%, 10%, 5%, and 3% contimination. 
To determine the g — r color offsets expected at constant 
J — H color between a sample of 997 dwarf stars (the 
number of Kepler EC hosts), we create a sample of 997 
dwarf stars by randomly selecting 997 stars from the orig- 
inal, not-intentionally contaminated control sample. As 
before, we compute the mean g — r color in bins of con- 
stant J — H color for each of the five samples described 
above. We plot the result of this calculation in Figure [BJ 

We find that between 10% and 30% of the stars in our 
control sample of stars with KIC-based log g > 4 would 
need to be giants to reproduce the 0.08 mag offset in g — r 
at J — H = 0.62 we observe between the hosts of Ke- 
pler EC and our original not-intentionally contaminated 
control sample. The level of contamination is at least 
a factor of a few lar ger than th e level of c ontamination 
found in the KIC bv lBasri et al.1 (|2011l ) and lBrown et al.1 
(|2011l) . Moreover, Figure [5] indicates that a giant star 
contaminated control sample would have redder g — r 
color at J — H = 0.54 than control sample of dwarf 
stars (very much like the likely hosts of Kepler ECs). 
We do not observe this effect in Figure [2] Indeed, the 
hosts of Kepler ECs are redder than the control sample 
at J — H = 0.54. At the same time, the presence of 
giant stars has no effect on the g — r color of our sam- 
ple at J - H = 0.22, 0.30, or 0.38. For that reason, the 
presence of giant stars in our control sample cannot ex- 
plain the significant color offsets we observe between the 
hosts of giant Kepler ECs and our original control sam- 
ple. We therefore conclude that contamination of our 
control sample by giant stars is unlikely to produce the 
color offsets we observe in Figure [2] and Table [TJ 

2.3. Metallicity Offsets from Color Offsets 

We use two methods to transform the g — r color off- 
sets we observe in Figure [5] and Table [JJ into approxi- 
mate metallicity offsets. First, we use dereddened fidu- 
cial sequences for the open clusters M 67 ([Fe/H] = 0.0 
and age » 4 Gyr) and NGC 6791 ([Fe/H] = 0.4 and 



age ~ 10 Gyr) in the Sloan photometric system from 
lAn et all (|2008[ ) plotted in Figure In this case, the 
fiducial sequences show that at constant 6 < M r < 7, 
the metal-enriched NGC 6791 fiducial sequence is w 0.1 
mag redder in g — r than the solar metallicity M 67 fidu- 
cial sequence. Though these two clusters have different 
ages, the stars with 6 < M r < 7 have main-sequence 
lifetimes much longer than 10 Gyr and therefore should 
still be on the main sequence in both clusters. Conse- 
quently, the « 0.1 mag offset in g — r between the two 
fiducial sequences is likely a result of the « 0.4 dex offset 
in [Fe/H] between the two clusters. Accordingly, at con- 
stant 6 < M r < 7 it seems that A [Fe/H] = 0.1 results in 
A(g - r) » 0.025. 

To transform the J — H colors of the points in Figure [2] 
and Tabled] into M r , we use a 2 Gyr solar metallicity 
Padova isochrone. To more precisely estimate the aver- 
age g - r offset between the M 67 and NGC 6791 fiducial 
sequences, we calculate the difference between the fidu- 
cial sequences at three points (M r — 5.9, 6.4, and 6.8 
corresponding to J — H = 0.46, 0.54, and 0.62) and av- 
erage the result. In this way, we find that in the interval 
6 < M r < 7, an offset of 0.1 dex in metallicity results 
in an offset of 0.035 mag in g — r. This method is only 
applicable to stars with 6 < M r < 7, so in Table [1] we 
only give the metallicity offsets indicated by this method 
for the three reddest J — H bins in our analysis. 

In order to transform a g — r color offset into a 
metallicity offset for more massive stars, we use Padova 
isochrones. In particular, we use two Gyr isochrones for 
two metallicities ([Fc/H] = -0.1 and 0.1) to determine the 
change in g — r color at constant J — H attributable to a 
0.2 dex increase in metallicity. We use the two isochrones 
to determine the offset in g — r at constant J — H that 
can be attributed to metallicity, then use that g — r off- 
set to convert the observed g — r offsets in Table [TJ into 
metallicity offsets (extrapolating if necessary). 

Interestingly, even in the presence of possible system- 
atics in our analysis (e.g., differential reddening, age dif- 
ferences, or the possible presence of giant stars in the 
control sample), it is reassuring that we recover the fact 
that solar-type stars that host giant planets are metal 
enriched relative to the field population. In Figure [5J we 
plot T c g and metallicity for both a volume-limited sample 
of stars (d < 20 pc) from the Geneva - Copenhagen Survey 
(GCS - iHolmberg et~aHl2007l [20091) an d confirmed giant 
planet hosts (e.g.. iWright et al.ll2011[ ). The horizontal 
lines give the average metallicities of the two popula- 
tions. The vertical arrows are at the T e g corresponding 
to J — H = 0.22 and 0.30. The arrows are anchored at 
the mean metallicity of the GCS sample, and their length 
corresponds to the metallicity offset we need to explain 
the differences in g— r color we observe in the giant planet 
bins at those J — H colors reported in Table [TJ Reassur- 
ingly, the metallicity offsets we compute based on g — r 
color offsets are in quantitative agreement with previous 
observations of metal-enhancement in giant planet hosts. 
In other words, our analysis is quantitatively correct in 
the regime where it is possible to compare with previ- 
ous results. That is suggestive that our result for the 
low-mass stars that host small-radius Kepler ECs is also 
reliable. 

3. DISCUSSION 
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We identified three properties of the Kepler ECs in 
Section 2: 

1. The population of solar- type stars that host giant 
ECs is 4-cr redder by 0.04 mag in g — r at constant 
J — H = 0.3 than a control sample of stars observed 
by Kepler but with no detected planets. 

2. The population of solar- type stars that host small 
ECs is similar in g — r at constant J — H = 0.3 to 
a control sample of stars observed by Kepler but 
with no detected planets. 

3. The population of low- mass stars that host small 
ECs is 4-cr redder by 0.08 mag in g — r at con- 
stant J — H = 0.62 than a control sample of stars 
observed by Kepler but with no detected planets. 

The first observation is consistent with the known 
correlation between solar-type host star m etallicity and 
likeli h ood of hosting a giant planet (e.g., iSantos et al.l 
12004 iFischer fc Valentil 120051) . The second observa- 
tion supports the tentative assertion that for solar- 
type stars in the solar neighborhood, metallicity is only 
weakly correlated with the likelihood of hosting a low- 
mass planet (e^ lUdrv et all 120061: iSousa et al.l 120081: 
iBouchv et al.ll2009l ). If the g — r color offset is the re- 
sult of a 0.2 dex enhancement in [Fe/H] of the low- mass 
hosts of small-radius Kepler ECs, then the third ob- 
servation supports the t entati ve assertion put forth in 
Schlaufman & Laughlin (2010) that M dwarfs that host 
low-mass planets are more likely to be metal rich than 
expected based on random sampling of the field M dwarf 
population. 

The observation that there is a correlation between 
host star metallicity and the presence of low- mass planets 
around low-mass stars is expected in the core-accretion 
model of planet formation. In that scenario, the cores of 
Neptune-mass planets (and possibly super-Earths) form 
and grow to several Earth masses before their parent 
protoplanetary disk is fully dissipated. To first order, 
the mass of a protoplanetary disk scales as Mdisk °c M* 
and the fraction of the disk in solid material scales 
as /solid oc Z*. The total amount of solids in a disk 
able to form the core of a Neptune-mass planet is then 
Msoiid /soiid-^disk oc Z*M*. The protoplanetary disk 
in which the solar syst em formed like ly had a solid mass 
Msoiid ~ 100 Mq (e.g.. lLissauerl ll993). If the probability 
of forming a planet is proportional to the total amount of 
solids in its parent protoplanetary disk, a late K dwarf 
with M* = 0.7 Mq would need to have a metallicity 
of [Fe/H] = 0.15 to have the same chance of forming a 
planetary system as a solar metallicity solar-type star. 
As a result, it may only be the metal-rich late K dwarfs 
that have protoplanetary disks with enough solid mass to 
form the several Earth-mass core of an ice giant or super- 
Earth. Though comparing stars at constant J — H color 
forces us to compare stars of slightly different masses, the 
magnitude of this effect on the planet formation process 
is likely an order of magnitude smaller than the effect of 
a 0.2 dex enhancement in [Fe/H]. 

The correlation between the red optical color of late K 
dwarfs and probability of hosting a small EC indicates 
that Kepler might boost its yield of planets by shading 
their late K dwarf sample to redder optical colors. In 



particular, focusing on stars in the range 0.5 < J — H < 
0.7 above the line 

g-r = 1.5[(J-H)-0.5] + l (2) 

might produce a larger sample of small ECs than a search 
that does not use red g — r color in its target selection. 

4. CONCLUSION 

We find that at J — H = 0.62, low-mass stellar hosts 
of small-radius Kepler ECs are 4-cr redder in g — r than 
a control sample of stars with the same J — H color but 
with no detected planets. This result is unlikely to be 
an artifact of reddening, age differences between the two 
populations, or the presence of giant stars in the con- 
trol sample. Stellar models and fiducial sequences in the 
Sloan photometric system of the Galactic open clusters 
M 67 and NGC 6791 indicate that a 0.2 dex enrichment in 
[Fe/H] of the EC host population can produce the offset. 
We suggest that small planets are preferentially found 
around metal-rich low-mass st ars (confirming the ten- 
tativ e correlation advocated in Schlauf man fc Laughlin! 
120101) . We also confirm that giant Kepler ECs are pref- 
erentially found around metal-rich solar-type stars and 
that the presence of small Kepler ECs orbiting solar-type 
stars does not seem to be strongly correlated with host 
metallicity. 
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Figure 1. Kepler exoplanet candidate (EC) host stars in a J — H — g — r color-color plot. We plot in red stars that host at least one 
giant EC (i.e., R p > 5 -R®), while we plot in blue stars that host an EC system with no giant planets. We plot in gray a control sample 
of 10,000 stars randomly selected from the ~150, 000 sta rs observed in Ql and Q2 of the Keple r mission that have no detected ECs. The 
black curve is a 2 Gyr, solar me tallicity Padova isochrone (Marigo_et al. 2008; Girardi ct al. 2010), and the orange curve is the M 67 fiducial 
sequence from|Xn ct al. (2003). We indicate with yellow arrows the affect of increasing metallicity 0.2 dex in [Fe/H] on g — r at constant 
J — H. Note though that the yellow arrows do not connect stars of constant mass, as a metal-enriched star will be about 5% more massive 
than a solar metallicity star at constant J — H color. We give approximate stellar mass and spectral type as a function of J — H color at 
the top of the plot. At J — H > 0.6 (typical of late K dwarfs), giant ECs become very rare relative to smalle r ECs. In other words, the 
Kepler ECs confirm the correlation between host stellar mass and frequency of giant planet occurrence (e.g., Johnson ct al. 2007, 120101 ; 
Borucki ct al. 2011b). Moreover, at J — H sa 0.62 characteristic of late K dwarfs, the effect of metallicity moves metal-rich stars to redder 
g — r at constant J — H. At that color, the population of K dwarfs that host ECs has a redder g — r color than the control sample (see 
Section 2.3). That offset is unlikely to be explained by differential reddening between the small EC host sample and the control sample, 
differences in age on the main sequence, or the presence of giant stars in the control sample (see Section 2.2). As a result, the sample of 
low-mass stars that host small ECs is likely more metal rich than the control sample of stars in the same range of J — H . 
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Figure 2. Left: mean g — r color as a function of J — H. Right: median g — r color as a function of J — H . We plot in red the result 
for the hosts of giant ECs, in blue the result for hosts of systems without a giant EC, and in gray the result for the control sample. The 
error bar on each point gives the 1-tr error on our estimate of the mean or median. We give approximate stellar mass and spectral type 
as a function of J — H color at the top of the plot. We give the number of control stars, giant EC hosts, and small EC hosts in each 
0.16 mag bin in J — H color at the bottom of the plot. For solar-type stars, the hosts of giant ECs have significantly redder g — r at 
constant J — H 0.3 (and therefore likely higher metallicity) than the control sample, reprodu cing the observed correlation between giant 
planet occurrence and host stellar metallicity (e.g., Santos et al. 2004; Fischer & Valcnti 2005). On the other hand, there is no significant 
difference between the average g — r color of small EC hosts and the control sample, confirming obser vations that for so l ar-type stars there 
is little correlation between stellar metallicity and the likelihood of hosting a low-mass planet (e.g.. IUdrv"et~a l. 2006; Sousa ct al.l l2008l ; 
IBouchv ct al. 2009|). At J — H ss 0.62 however, the hosts of small ECs are significantly redder (and therefore likely more metal rich) than 
the stars in the control sample with no observed ECs. We report the numerical values plotted above in Table [T] This result is a natural 
consequence of t he core-accretion model of plan et formation and confirms at much higher statistical significance the tentative hint of that 
relation noted in Schlaufman &z Laughlin (2010). 
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Figure 3. Color-color and color-magnitude diagrams showing the distribution of stars in the control sample in gray, stars that giant 
Kepler ECs in red, and stars that host small-radius Kepler ECs in blue. In all cases, the control sample traces the same parameter space 
as the sample of EC hosts. The EC hosts that are outliers in each color-color plot may be unequal mass binary star systems or systems 
with poorly-estimated reddening. 
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Figure 4. Left: estimated color excess E (B — V) of Kepler EC host stars from the KIC as a function of apparent Kepler magnitude K p . 
Right: equatorial coordinates of the stars used in this analysis. In both panels, we plot in red stars that host at least one giant EC, while 
we plot in blue stars that host an EC system with no giant planets. We plot in gray the control sample. The average color excess of both 
the control sample and the giant EC host sample is 0.12, while the sample of EC hosts with no giant planet has an average color excess of 
0.10. The sharp edge in E (B — V) at K p = 14 is likely a selection effect: apparently bright solar-type dwarf stars are relatively rare, so 
they were included in the ~150, 000 stars to be searched for transiting planets in Ql and Q2 of the K epler mission regardless of reddening. 
On the other hand, faint solar-type dwarfs stars are numerous, so only the least reddened stars were included in the sample to be searched 
for transiting planets. The spatial distribution of EC hosts on the sky is similar to the distribution of stars in the control sample. In short, 
the EC host sample is subject to the same reddening as the control sample. Therefore, any optical color offsets between the control sample 
and the EC host sample are unlikely to result from differential reddening. 
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Figure 5. Timescale for tidal migration to move a planet on an initially circ ular orbit to one stellar r adius, a reasonable approximation 
to the timescale to tidal disruption. In this case, we assume that Q' t = 10 6 , aslSchlaufman ct al. (2010) found that the average Q'^ in the 
Kepler EC host sample was 10 6 < Q'„ < 10 7 . The median period of giant ECs is 13 days, while the median period of small-radius ECs is 
10 days. As a result, tidal evolution does not significantly affect the bulk of the Kepler EC, so there is no reason to expect the hosts of 
Kepler ECs to be preferentially young stars. 
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Figure 6. Effect of giant star contamination in the control sample. We find that between 10% and 30% of our control sample with 
KIC-based log g > 4 would have to be misclassificd giants to explain the the 0.08 mag g — r color offset we observe at J — H = 0.62 between 
the control sample and the hosts of small-radius Kepler ECs. This level of contamination is at least a factor of a few larger than the level 
of contamination expected in a sample of stars with KIC-based logg > 4 (Basri ct al. 2011; Brown ct al. 2011). Contamination by giant 
stars cannot explain the significant color offset we observe at J — H = 0.30 between the control sample and the hosts of giant Kepler ECs. 
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Figure 7. Dereddened fiducial sequences for the op en cluster s M 67 ([Fe/H] = 0.0 and age m 4 Gyr) and NGC 6791 ([Fe/H] = 0.4 and 
age ~ 10 Gyr) in the Sloan photometric system from An ct al. (2008]). At M r si 6.8 (corresponding to J — H = 0.62 characteristic of late 
K dwarfs), the NGC 6791 fiducial sequence is 0.1 mag redder in g — r than the M 67 fiducial sequence. This offset is unlikely to be due to 
the age difference between the two clusters, as the main sequence lifetime of a late K dwarf is much longer than 10 Gyr. In other words, 
at J — H = 0.62, every 0.1 dex increase in [Fe/H] corresponds to an increase in g — r of about 0.025 mag. 
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Figure 8. Scatterplot of ef fective temperature T e gf and metallicity for both a volume -limited (d < 2 pc) sample of stars from the Geneva- 
Copenhagen Survey (GCS - Holmbcrg ct al. 2007, 2009) and giant planet hosts from Wrieht ct al. (201f). We plot in gray the GCS stars 
and in red the hosts of giant planets; the horizontal lines give the average metallicities of the two populations. The yellow vertical arrows 
give the metallicity offsets suggested by the g — r color offsets we observe between the Kepler giant EC hosts and the control sample at 
J - H = 0.22 and 0.30 reported in Table [T] Indeed, though our analysis is possibly imprecise and affected by several systematic effects 
(e.g., differential reddening, age differences, or the possible presence of giant stars in the control sample), we reproduce both qualitatively 
and quantitatively the known metallicity offset between solar neighborhood stars that host giant planets and those that do not host giant 
planets. 
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Table 1 

Optical Color Offsets and Equivalent Metallicity Offsets 



J -H 

(mag) 


(Mq) 


A( S - r) 
(mag) 


Significance 


A[Fe/H] a 
(dex) 


A[Fe/H] 
(dex) 


b A( S - r) 
(mag) 


Significance 


A[Fe/H] a 
(dex) 


A[Fe/H] b 
(dex) 


0.22 


1.23 


0.01 


1.70 




0.11 


0.03 


2.06 




0.22 


0.30 


1.05 


0.01 


1.14 




0.05 


0.04 


4.44 




0.22 


0.38 


0.93 


0.01 


1.91 




0.07 


0.02 


2.59 




0.14 


0.46 


0.85 


0.06 


1.72 


0.17 


0.39 


0.08 


4.11 


0.23 


0.54 


0.54 


0.78 


0.02 


1.03 


0.06 


0.09 


0.02 


0.69 


0.06 


0.09 


0.62 


0.73 


0.08 


4.00 


0.24 












Note. 


— Columns three throuj 


;h six correspond to the stellar hosts of EC 


systems with no 


giant ECs, while columns seven 


through ten 



correspond to the stellar hosts of EC systems with at least one giant EC. 
a Metallicity offset based on M 67 and NGC 6791 fiducial sequences from lAn et~aLl pOOST ). They are only applicable for stars with 
J — H > 0.46. 

b Metallicity offset based on Padova isochroncs (Marigo et al. 2008; Girardi ct al. 2010). They are only applicable for stars with J—H < 0.54. 



